Nonlinear screening of charges induced in graphene by metal contacts 
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To understand the band bending caused by metal contacts, we study the potential and charge 
density induced in graphene in response to contact with a metal strip. We find that the screening 
is weak by comparison with a normal metal as a consequence of the ultra-relativistic nature of the 
electron spectrum near the Fermi energy. The induced potential decays with the distance from the 
metal contact as x _1//2 and x _1 for undoped and doped graphene, respectively, breaking its spatial 
homogeneity. In the contact region the metal contact can give rise to the formation of a p-p' , n-rt , 
p-n junction (or with additional gating or impurity doping, even a p-n-p' junction) that contributes 
to the overall resistance of the graphene sample, destroying its electron-hole symmetry. Using the 
work functions of metal-covered graphene recently calculated by Khomyakov et al. [Phys. Rev. B 
79, 195425 (2009)] we predict the boundary potential and junction type for different metal contacts. 
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I. INTRODUCTION 

Graphene's unique electronic properties arise from the 
ultra-relativistic character of the electron spectrum near 
the Fermi energy that leads to many unusual physical 
effects The high mobility and low electron density 
(compared to normal metals) of this two-dimensional 
form of carbon suggest various possibilities for using 
single and multiple graphene sheets to make electronic 
devices pffio] p or example, one can use external gates to lo- 
cally change the doping of graphene and so design p-n or 
p-n-p' j unctions How charge inhomogeneities are then 
induced in graphene by the gate electrodes or by charged 
impurities has been studied theoretically in Refs. [T2HT51 

Recently, we have shown how a graph ene s heet ad- 
sorbed on a metal is charged by the metaP^^ suggest- 
ing an alternative way of making p-n junctions by putting 
weakly-interacting metal strips on graphene. The (out- 
of-plane) charge transfer between graphene and the metal 
is determined by (i) the difference between the work func- 
tion of graphene and the metal surface, and (ii) the metal- 
graphene chemical interaction that creates an interface 
dipole which lowers the metal work function. Depositing 
a metal strip (electrode) of finite width on the graphene 
sheet will additionally result in an in-plane charge trans- 
fer from the graphene region covered by the metal to 
the "free" graphene supported by a dielectric, driven by 
the difference between the work functions of the metal, 
Wm, of metal-covered graphene, W, and of free-standing 
graphene, Wq. 17 An electrostatic potential will then be 
induced across the graphene sheet, leading to band bend- 
ing and the formation of a p-p' n-n' or p-n junction at 
the contact area as illustrated in Fig. jipilSI 

In this paper we study how graphene screens charges 
transferred from the metal- graphene contact to "free" 
graphene, creating an electrostatic potential barrier at 
the contact region, see Fig. [T] This problem is of fun- 
damental importance since the physics of contacts can 
have a significant effect on the transport properties of an 



electronic device, in particular, graphene-based devices 
with "invasive" electrodes, as recently demonstrated in 
Ref. 18. A widely-used picture of a metal-graphene 
contact often assumes a quite sharp potential step at 
the contact.^ However, we find that the screening in 
graphene is strongly suppressed leading to a large-scale 
inhomogeneity of the electrostatic potential across the 
graphene sample. The contact effects can also result in 
charge inhomogeneities of different polarities (an abrupt 
p-n junction) at the contact region. These abrupt p- 
n junctions are qualitatively different from the gate- 
induced p-n junctions recently reported in Refs. I6|ll|14l 

The contact phenomenon studied in the present work 
stems from long-range electric fields originating from 
charge transferred between a graphene sheet and a metal 
electrode. Recently, Barraza et al. have obtained a con- 
tact potential from density functional theory calculations 
of graphene ribbons with a length of < 13.6 nm contacted 
with aluminium electrodes. 20 To calculate the contact po- 
tential profile one applies periodic boundary conditions, 
i.e. periodic arrays of metal-graphene junctions. Using 
such results to describe a single metal-graphene junction, 
one assumes that the periodic images do not interact, or, 
in other words, that the contact potential profile is con- 
fined to a narrow region around the junction. Indeed 
the phenomenological model proposed in Ref. [20] sug- 
gests an exponential-like decay of the potential induced 
in graphene by a metal contact. Our results show that 
this does not apply to a single metal-graphene contact, 
where the induced potential is in fact long range. 

In Sec. [IT] we describe the Thomas- Fermi approach used 
to study the band banding caused by metal contacts in 
graphene. In Sec. |III| we derive the induced potential and 
charge density for undoped, doped and gated graphene, 
and make material-specific predictions of the boundary 
potential and junction type for different metal contacts. 



Our conclusions are presented in Sec. IV 
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II. THOMAS-FERMI APPROACH 



potential induced in graphene (x > 0) is 



We model a graphene- metal contact as shown in Fig.[T] 
Metal contacts are often applied as a bilayer; a (thin) 
layer of metal 1 takes care of a good adhesion to 
graphene, and a (thick) layer of metal 2 forms the elec- 
trode. We allow for electrostatic doping by means of 
a back gate electrode, and model chemical doping by a 
uniform background charge a. The resulting screening 
potential in graphene is calculated using density func- 
tional theory (DFT) within the Thomas-Fermi (TF) ap- 
proximation. The latter has proven to have a wide range 
of validity for describing screening in graphene i n 
TF theory the charge density cr(x) induced in graphene 
for x > is given in terms of the local chemical poten- 
tial ji(x) as <j{x) — —eDojj J (x)\jj J (x)\/2A^ where Do = 
0.09/(eV 2 unit cell), A = 5.18 A 2 is the unit cell area of 
graphene, and e > the elementary charge.^ The elec- 
trostatic potential V(x,z) and the local chemical poten- 
tial are related as ji(x) = /j J p — V(x) with V(x) = V(x, H) 
and the graphene sheet is located at z = H, where H is 
the distance between the graphene sheet and the back- 
gate electrode, see Fig. [T] The chemical potential of the 
entire system, /ip, is set to zero for undoped graphene. 

Using the Poisson equation and translational symme- 
try in the y direction, (i.e. assuming an infinitely wide 
graphene sample), the TF equation for the electrostatic 




FIG. 1: (Color online) Schematic picture of a metal strip 
(electrode) on graphene. V(x,z) is the metal-induced electro- 
static potential, where V(0,z) is fixed by the work functions 
Wmi, Wm2, Wg and W, of the bottom and top metal layers, 
graphene and graphene-covered metal; hi, and H are the 
thicknesses of the bottom and top metal electrodes and the 
insulating substrate, and Lq and L are the lateral dimensions 
of the electrode and the graphene sample. <J> g is the back- 
gate voltage; Ko and ki are the permittivities of the dielectric 
media above and below the graphene sheet. A typical exper- 
imental setup has Lq < L ~ 1 /mi, H ~ 300 nm < L, and 
hi < h 2 < H. 



^V(x,z) = e[a{x) + ° ] S(z-H), 
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where a(x) = — eA[/ip — V(x)]\fi^ — V(x)\ with A = 
Do/2 A. The effects of the substrate, the lattice potential 
and of the filled band electrons are taken into account 
via an effective dielectric constant, ft, and a compensat- 
ing charge a = eA /ip |mf I ! 1S | 21 1 The boundary conditions 
for V(x, z) are imposed by the potential at the contact 
V(0, z) = Vc(z) and at the gate V(x, 0) = V g (x). Assum- 
ing graphene is a nearly perfect metal, i.e. V(x) « const, 
which will be justified below, the boundary conditions 
are given by^HH 
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where z > and 7r/2</3<7ris the angle between the 
"free" graphene sheet and face of the metal electrode; 
W^mi? Wm2, Wq and W are the work functions of the 
bottom and top metal layers, graphene and graphene- 
covered metal, respectively. <l> g is the backgate voltage. 

Using a Green function approach, Eq. Jl]) can be for- 
mulated as a nonlinear integral equatiorP^T 
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Using Eqs. ([2| and (|3| the boundary potential term can 
be written as 
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where the boundary potential constants Vbi and Vb2 are 

Vm = ±(W-W G ), V B2 = ^(W M -W G ), (7) 

and W and W& are the work functions of the graphene- 
covered metal (Ml) surface as calculated in Ref. 17 using 
the DFT within the local density approximation, and 
that of free-standing graphene, respectively. The influ- 
ence of the electrode depends upon the geometry and 
the electrode work function Wm — Wmi for /12 = and 
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Wm = Wm2 for I12 ^> hi. The parameter f3 of the con- 
tact geometry is tt/2 for x <C h = hi + h<i and tt for 
x ^> h. At a distance x ~ h the two solutions of Eq. Q 
obtained with these two values of f3 serve as lower and 
upper bounds for the screening potential V(x). 

In terms of the sine Fourier transform f(k) = 
\/2/it J °° dx f(x) sin(fcx), Eqs. Q and (J^J) reduce to 



Ve(fc,* g ) 
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where a(k) - 
(2/nk 2 ) 1 / 2 a. 



' dx a(x) s'm(kx) and a(k) 



III. RESULTS 



We have been unable to solve Eq. Q analytically, but 
we have obtained a numerical solution, as well as ac- 
curate, approximate, analytical solutions. We start by 
studying the solution at a large distance from the metal 
contact assuming that V(x) — » Vdo = const for x — » oo. 
The classical limit of a perfect metal, i.e. V(x) = Vqq for 
all x, implies V(k) = {2/nk 2 ) 1 ^ 2 V 00 in Eq. ([8]), allowing 
the asymptotic charge density to be written as 



e [cr(x) 
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which is derived from Eqs. ([8| and ([9| with the use 
of the inverse sine Fourier transformation. This result 
is valid for any two-dimensional metallic system, in- 
cluding multiple layers of graphenePS However, to ob- 
tain the screening potential V(x) from the charge den- 
sity a(x) one has to use the appropriate density of 
states. For instance, for undoped graphene, /ip = 
(implying a = 0) and the asymptotic potential is 
Voo = sign(— e^> g )(£o^|^ ) g|/eAi^) 1 / 2 . For a single layer 
of graphene the screening potential can be expressed in 
terms of the charge density in the general form 



V(x) = fi F + sign(a) 



1*0*01 

eX ' 



(11) 



where the charge density is given by Eq. (10) for given 
boundary potential constants (Vbi, Vb2)> gate voltage 
(3>g), substrate thickness (H) and doping level (//f)- 

In the ungated limit (H — > oo), the screening potential 
given by Eqs. (10) and (11) becomes 



V(x) = fi F + sign(a) 1 



V^b|Vb| 



x/l a 



(12) 



1.0* 
0.8 ^ 



^ 0.6 



0.4 

0.2 1 
0.0 : 





1.0: 




U.j : 


> 




CD 


0.0: 


b 


-0.5 : 




-1.0 




A pn 



Ti/Au 



20 40 60 80 100 







20 



40 60 



80 



100 



FIG. 2: The screening potential for undoped graphene: 
num eric ally exact solution (solid line), variational solution 
Eq. ^ (bold dots), V(x) = V B (x/l s )~ 1/2 (dotted line) and 
V{x) — Vb (1 + x// s ) -1 / 2 (dashed line). Inset: the charge 
density a(x) = eXV(x)\V(x)\, Eq. |l0|, for the Ti/Au elec- 
trode for <J> g = — 15 V (dashed line), V (solid line) and 15 



V (dotted line); Vb = 0.19 eV, 
A)" 2 , H 300 nm, l s ~ 1 nm. 



2.5, A = 8.69 x 10" 3 (eV 



where l s = hv/7ra\VB\ is a scaling length, Vb = Vbi + 
^B2, hv = ±e k /\k\ = (ttX)-^ 2 = 6.05 eVA, 16 and a = 
e 2 / AitEonhv = 2.38/ft the "fine-structure" constant in 
graphene. 



A. Undoped graphene 

For undoped graphene (/ip = 0) the induced potential 
decays asymptotically as V(x) ~ x -1 / 2 , i.e. the screening 
is strongly suppressed as compared to a normal 2D metal, 
where V(x) ~ x~ x . The charge density a(x) ~ V 2 (x) ~ 
x~ x behaves, however, as in a 2D metal.^ 

Numerical solution of Eq. Q shows that Eq. ( [l2| ) is 
quite accurate for x ^> l s . A simple interpolation ex- 
pression for the screening potential of undoped graphene 
(/xf = 0) can be obtained by replacing x/l s with l + x/l 8 
in Eq. ( 12 ). The screening potential then satisfies the cor- 

and at x 



rect boundary conditions both at x = U and at x = oo, 
and is close to the exact solution also for intermediate x, 
see Fig. [2p^ The best variational solution we found is 



V{x) 



Vb 



(^x/l s + p 2 ft) V2 (x/l s + /3r 2 ) 



\V4 



(13) 

with f3i = 0.915 and f3 2 = 0.128. The difference between 
the numerical and variational solutions is < 1.4% for x < 
0.1 l s and < 0.2% for x > 0.1 l B . 

From the general criterion for validity of the semiclas- 
sical approximation, \d\(x)/dx\/2ir <C 1, where X(x) is 
the de Broglie wavelength, 26 one can derive a condition 
defining the range of validity of the potential obtained 
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TABLE I: Vb, Eq. ([7|, is the boundary potential for double layer electrodes with hi <C fe, and for two contact geometries with 
f3 = 7r/2 and 7r. Wm, Wg = 4.48 eV and W are the work functions calculated for close-packed surfaces of the clean metals, for 
free-standing and adsorbed graphene, respectively, see Ref. 1171 The doping sign of adsorbed and "free" graphene corresponds, 
respectively, to the sign of W — Wg and Vb = Vbi + Vb2- 
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from the TF model: \dV(x)/dx\ <C V 2 (x)/hv. The 
latter condition also implies that the in-plane electric 
field is small compared to the field perpendicular to the 
graphene sheet, E x <C E z (as E z ~ a ~ V 2 ), which com- 
pletes the proof that graphene behaves as a nearly perfect 
metal within the TF approximation. Using Eq. (13), and 
taking a typical boundary potential Vb ~ 0.5 eV (Table [i]) 
and effective dielectric constant k ~ (k,q + k>i)/2 ~ 2.5 
for graphene on a Si02 substratepH we find that the TF 
theory is valid for x ^> a with a the lattice parameter of 
graphene. Using high k substrates will weaken this con- 
dition, making the TF results valid up to the immediate 
proximity of the metal contact. 

Recently we described how depositing graphene on a 
metal surface leads to charge transfer between the metal 
and graphene ! 16 * 17 * The graphene device shown in Fig. [I] 
consists of two regions, metal-covered graphene [x < 0) 
and "free" graphene (x > 0). The sign and level of doping 
of metal-covered graphen e is fix ed by the first metal layer 
(Ml), AE F = W - W G I =7 The doping of graphene 
near the contact, a <C x <C H/tt, is determined by the 



boundary potential Vb, see Eq. ( |13| ), which depends on 
the work function of the graphene-covered metal (Ml) 
and the top metal layer (M2), see Eq. 0. If sign(V B ) ^ 
sign(A^p) an abrupt p-n junction sh ould form close to 
the contact at x ~ a, see inset in Fig. [2] 27 Vb and AE-p 
having the same sign but a different size results in p- 
p' or n-n' junctions. Such junctions break the electron- 
hole symmetry in graphene, giving rise to a n asymm etric 
resistance as a function of the gate voltage P 8 * 2 I 29 * 

At a distance x ^> H/tt the induced charge is deter- 
mined by the top metal layer ("Metal 2"), Vb2 ~ Wm2 — 
Wg, whereas far from the contact, x ^> max(Lo,/i), 
graphene is unaffected by the metal electrode. Using 
the work functions calculated in Ref. [171 the boundary 
potential and junction type expected for different metal 
contacts are listed in Table [TJ assuming the metal surfaces 
are clean. Because of the sensitivity of work functions 
to surface contamination, the contact potentials may be 
different for experiments that are not performed under 
ultrahigh vacuum conditions.^ 




FIG. 3: (Color online) The screening potential for doped 
graphene, Eq. |l2| : /xf/Vb = —0.12 (dot-dashed green/gray 
line) and —0.25 (dot-dashed black line); uf/Vb = 0.12 
(dashed green/gray line) and 0.25 (dashed black line). As 
a reference we also show the screening potential for undoped 
graphene, /xf = 0, (solid line). Inset: the charge density for 
doped graphene. 



B. Doped graphene 

Real graphene samples are often doped with some form 
of charged impurities. Eq. ( [l2| ) shows that the screen- 
ing potential is relatively unchanged by doping for a 
graphene region where |V(x)| ^> |/xf|, and asymptoti- 
cally behaves as V^|T^|Z s /2|/iF|^ f° r 1^0*01 ^ ImfI? see 
Fig. [3j i.e. the screening is enhanced by doping. 

Eq. ( [T2] ) reveals an interesting effect related to the 
sign of doping. The screening potential is different for 
graphene doped with electrons and holes due to a p-n 
junction formation at x pn = l s (^f/^b) -2 for Mf/^b > 
as shown in Fig. [3] The formation of the p-n junction 
for /iF /Vb > is caused by the competition between the 
charges of opposite sign induced by the metal contact 
and the charged impurities, respectively. The screen- 
ing charge and potential for x < x pn (x > x pn ) are 
constrained by the metal contact (the charged impuri- 
ties). The screening at x ~ x pn is then suppressed, which 
gives rise to increase of the screening potential at the p-n 
junction as compared to the case of u f /Vb < OP For 
I^y/Vb < the charges induced in graphene by the metal 
contact and the charged impurities are of the same sign so 
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the overall effect leads to enhancement of the screening, 
reducing the induced potential. The asymmetry reaches 
its maximum at a critical doping level |/ip| ^ 0.1 |Vb| 
and vanishes upon increasing (decreasing) /ip, i.e. for 
\u F \ > 0.1 |y B | Of| < 0.1 1 Vb I)- This effect should be 
observable in transport measurements where the impu- 
rity doping is gradually changed from n- to p-type!^ 



electrodes (Ni, Co, Pd, Ag, Cu, Ti/Pd and Ti/Au), see 
Table [I] and Fig. [2] This is possible because the sign 
of the graphene doping beneath the electrode, near the 
contact and far from the contact is independently fixed 
by the bottom metal layer, the contact potential and the 
gate voltage (impurity doping), respectively. 



C. Gated graphene 



IV. CONCLUSIONS 



One can also create p-n junctions by applying a gate 
voltage (Fig.[l)PE2EaEq. ([ToJ shows that the first ef- 
fect of gating graphene is to cause a constant shift of the 
induced charge density. Secondly, the charge depletion 
width can be varied by changing the gate voltage since 
Voo ~ l^gl 1 / 2 . This accounts for the electric field ef- 
fect on the work function of graphene.^ Thirdly, a gate 
voltage applied to graphene creates qualitatively differ- 
ent junctions (n-n', p-p' , p-n or n-p) at the contact re- 
gion depending on the sign of the voltage. This will then 
lead to an asymmetry in the transport characteristics of 
the graphene device rather similar to the contact-induced 
abrupt junctions discussed in the last two paragraphs of 
Sec. HITAl 

The point x pn that separates p and n-doped regions 
of graphene, is given by <j(x pn ) = 0. Eq. (10) yields 
itXpn/AH w (Vb — Voo)/e<& z . The screening potential 
near the transition point can be calculated using TF the- 
ory. The semiclassical description, however, breaks down 
exactly at there are no screening charges at 

this point. The quantum corrections to the TF the ory 
are, nevertheless, relatively small in graphene! 12 * 14 * 15 * 

According to our analysis, gating graphene (or impu- 
rity doping) can create an n-p-n' junction for some metal 



We have studied the electrostatic barrier formed in 
graphene in response to a metal strip in contact with 
the graphene sheet. By comparison with conventional 
metals, the screening in graphene is strongly suppressed: 
the induced electrostatic potential decays weakly with 
the distance from the metal contact as V(x) ~ x -1 / 2 
and ~ x~ x for undoped and doped graphene, respec- 
tively. This leads to a substantial space charge region 
in graphene, breaking its spatial homogeneity. The lat- 
ter has been recently observed by scanning photocurrent 
microscopy, 8 and mi ght al so be seen for graphene in the 
quantum Hall regime pilMJ ^he contact effects also result 
in the formation of a p-p' , n-n' or p-n junction at the 
near contact area that breaks the electron- hole symme- 
try and contributes to the contact resistance! 18 * 28 * 29 * 32 * 33 * 
We predict that n-p-n' junctions can be realized by gat- 
ing graphene or by impurity doping in combination with 
specific metal electrodes. 
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